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Abstract 

Nonequilibrium phenomena of the phase transitions are studied. It is 
shown that due to finite relaxation time of the particle distributions, the 
use of scalar background dependent distribution functions is inconsistent. 
This observation may change the picture of rapid processes during the elec- 
troweak phase transition, like subcritical bubble formation and propagation 
of bubble walls. 



*riotto@tsmil9. sissa.it. Address after November 95: Theoretical Astrophysics Group, 
NASA/Fermilab, Batavia, IL60510, USA. 

tvilja@utu.fi. 



The possibility that the Universe went through a series of phase transitions during 
its expansion and cooling down from temperatures close to the Planck scale has been 
intensively investigated in the past [T|]. It is a common hope that many of the current 
questions of cosmology can be answered by studying the nontrivial dynamics of the ap- 
proach to equilibrium in complex systems. Nevertheless, despite their immense relevance, 
only very recently more effort has been made to understand nonequilibrium phenomena 
occured in the early Universe. 

Thermalization, reheating and relaxation are nonequilibrium phenomena that play 
crucial role in the completion of the inflationary era, which is thought to solve the horizon 
and the homogeneity problems [fjj of modern cosmology. Recent investigations on the 
non-linear quantum dynamics of scalar fields have implications for the reheating and 
reveal that particle production appears to be significantly different from linear estimates 
due to the time evolution of the inflation field |J. The quantum non-linear effects lead to 
an extremely effective dissipational dynamics and particle production even in the simplest 
self-interacting theory where the single particle decay is kinematically forbidden. Also 
they might help in alleviating the so-called Polony M problem concerning flat directions 
for some of the moduli fields in string theories. 

At the electroweak scale, the focus has been in generation of the baryon asymmetry 
during a first-order phase transition where the SU(2)l ® U(1)y symmetry is broken 
down to U(l) em ||. It is currently believed that nonequilibrium conditions are a crucial 
ingredient for electroweak baryogenesis, even though there are certain questions related 
to the reliability of the perturbative expansion for weak transitions 0. Moreover, the 
mechanism of the weak first-order transitions @ is still an open question. Consider 
models with double-well potentials where the system starts localized on one minimum. 
For sufficiently weak transitions subcritical bubbles of the other phase could be thermally 
nucleated, giving rise to e.g. an effective phase mixing between the two available phases 
before cooling down to the tunneling temperature at which critical bubbles are expected 
to be nucleated. This may have dramatic consequences for any electroweak baryogenesis 
mechanism. 

So far, however, most done to study the nonequilibrium aspects of phase transitions, 
e.g. the influence of subcritical bubbles on a first-order electroweak phase transition, is 
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strictly related to the equilibrium finite temperature field theory || and effective potential 
which, by its very definition, is only adequate to describe equilibrium situations. Two 
equivalent methods are usually applied: the imaginary-time formalism (in Euclidean 
time) or the real-time formalism with equilibrium distribution functions. They are used 
to obtain the thermal part of the equilibrium effective potential by integrating out all 
the fermionic and bosonic degrees of freedom of the theory with c -dependent masses. 

In the equilibrium case the finite temperature part of the one-loop effective potential 
reads || (besides an additive c -independent constant) 

(2^p L ^AM> W 



where uj^rrii) = y |k| 2 + mf(<f) c ), mj(0 c ) is the eigenstate of the mass matrix correspond- 
ing to the i-th degrees of freedom in a c -dependent background. Note, that /* is 
the usual equilibrium Fermi-Dirac or Bose-Einstein distribution function. The effective 
potential V eq (T, <p c ) represents the equilibrium free energy density as a function of the 
classical order parameter <fr c and can be used to determine the nature of the phase transi- 
tion and static quantities like the critical temperature, but it is not the appropriate tool 
for the description of real-time processes, which are crucial to understand the mechanism 
by which the system approaches the equilibrium . 

The aim of the present Letter is to point out that the use of the equilibrium finite 
temperature effective potential becomes completely unreliable when the thermalization 
time r t h of the degrees of freedom getting a c -dependent mass in the thermal bath is 
larger than the typical time scale C /0 C , I.e. when the macroscopic order parameter 
4> c does not vary slowly in comparison with r t h. More generally, we can state that the 
equilibrium effective potential is not trustable at times t smaller than the thermalization 
time T t h in a c -background changing in space and time at a generic instant t = 0. 

To avoid any confusion, we observe here that what we call thermalization time should 
be understood as the time needed for the degrees of freedom with a c -dependent mass 
to feel the change of the background <p c and to relax their distribution functions to their 
equilibrium values /* q with the appropriate mass mj(0 c ) [|ll|]. Thus, for instance, the 
thermalization time for fermions is dictated by the Yukawa interactions with the scalar 
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field and not by the strong interactions with gluons in the thermal bath, which are 
much faster. 

In this paper we propose a nonequilibrium approach to describe the properties of the 
system in the specific limits depicted above and apply it to the Standard Model (SM) 
effective potential. We will also briefly discuss some physical situations in which our 
simple, but crucial, statement might have interesting consequences. 

Let us first briefly clarify with a simple example our observation. Suppose a universe 
filled only with a real scalar field and a fermion field ip whose Lagrangian density reads 

£ = \ (d^ f + i$W>~ V(4>) - g<m. (2) 
The zero temperature tree level potential V{4>) is given by 

v(4>) = -\W + \\^. (3) 

At high temperatures the effective potential for the classical field (0) = C has the form 

V(TA c ) = ±m\T)<j ) t+ 1 -\(T)<pj, (4) 

where 

m\T) = V + (j + y) T\ (5) 

is the plasma mass for the <f> c - field whereas A(T) receives only logarithmic one-loop cor- 
rections. Above the critical temperature T c , the free energy of the system is minimized 
by 0c = and the Fermi-Dirac and Bose-Einstein distribution functions for the fermionic 
and bosonic degrees of freedom of the heat bath are the equilibrium ones 

/; q (k,0 c = O) = (l + e^l) _1 

f;\k, C = o) = (Ww 2 w ( T) _ ^ - 1 (6) 

Let us now imagine that at a generic time t = in a certain region of space the 
vacuum expectation value becomes different from zero, C = c (x, t). This can happen 
because in that region the system suffers a thermal fluctuation or simply because the 
temperature of the thermal bath has cooled down to the critical temperature T c , C = 
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becomes an unstable point and C starts to roll down. The equation of motion for c (x, t) 
reads 

□0 C x, t + — * = 0, 7 

0<Pc 

where V((p c ) depends upon the distribution functions and through the one-loop 
corrections. The question is now whether we can use, for instance, the equilibrium dis- 
tribution function f^((p c ) with a c -dependent mass m^,(0 c ) = g(f) c to calculate the 
contribution of the fermionic degrees of freedom to the one-loop effective potential. 
To answer the question one should solve the the Boltzmann equation 

d t u + ±-vu + k-d k u = c[u\, (8) 

where C [/,/,] is the collision operator of the fermions. We define r t h is the thermalization 
time for the fermionic degrees of freedom, i.e. the time needed for the fermions of the 
thermal bath to response to the change in the background <ft c . By dividing for small 
times < 5t < r th the distribution function into two parts = /^ q (0 c = 0) + 5f^ the 
lowest order correction Sf^ is determined by the collision operator with the distributions 
replaced with /^ q (0 c = 0). The collision operator can be written in the form 

(2n) 3 E f (q) J^fE b {k) ] 
The dispersion relations are determined by rapid forward scatterings and thus the energy 



conservation relation is now given by y |k| 2 + m? + 3A0 2 = J |p| 2 + g 2 (f) 2 + y |q| 2 + g 2 <\> 2 - 
Applying it, it is straightforward to convince oneself that the square bracket term in 
the collision operator given by (||) gives a contribution of the order g 2 or A. Therefore, 
because also the matrix element is of the order g 2 , one obtains the result that 

ifit^r 1+0 ^) +0 ^)- (10) 

On the contrary, expanding the c -dependent equilibrium distribution 



around (fi c = 0, one should get 



fMc) = ( 1 + e^vW+^V 1 (11) 
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Tft^r l+ °^- (12) 

It is therefore clear that for times smaller than r t h, the distribution function is far 
from being equal to f^(4> c ) and, as a consequence, the use of the equilibrium one-loop 
effective potential is not appropriate to describe the dynamics of the system. 

This example contains the simple, but essential feature that thermalization takes 
time and, therefore, the degrees of freedom of the thermal bath can not always follow the 
change of the background <fi c . Indeed, thermalization requires real scattering processes 
and therefore is usually fairly slow. Forward scatterings, instead, do not change the 
distribution functions of particles traversing a gas of quanta, but modify their free dis- 
persion relations. This remains true also in the case of a nonequilibrium system. Forward 
scattering manifests itself, for example, as ensemble and scalar background corrections 
to the particle masses. Since the forward scattering rate is usually larger than the non- 
forward rate, non-equilibrium ensemble and scalar background corrections are present 
even for times smaller than the thermalization time. Afterwards, non-forward reaction 
rates become active and the system thermalizes. 

Technically, the thermalization rate 7 t h = t^ 1 is related to the imaginary part of 
the two-point function via j t h — Im r' 2 '(w,k)/w. Here u = u;(k) is the solution to the 
dispersion relations where the scalar dependent mass is involved. 

To decide how fast each particle species do thermalize, one should calculate the imag- 
inary part of the two point function for each particle species separately. It is physically, 
however, clear that the rates are essentially the same than their contributions to the Higgs 
thermalization rate. Regardless what the actual rates are, the concept of non-equilibrium 
at short times remains. In the Standard Model, at one loop, the imaginary part of the 
Higgs propagator receives contributions from decay and absorption (emission) processes. 
Absorption (emission) is always proportional to the difference between the distribution 



functions of the two external particles in the final and initial state [12 



At high temperature limit, the thermalization rate of leptons with the Higgs back- 
ground are negligible because of their small Yukawa couplings. Also for the quarks the 
thermalization rates with the Higgs background are small due to small Yukawa couplings 
apart the top quark which has the rate of the order j t ~ 10 -2 T or smaller. At two 
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loops kinematical constraints no longer exist. The leading, purely bosonic, contribution 
has been estimated in ref. |T2|] from the gauge boson loops and the thermalization rate 



for the gauge bosons with the Higgs background turns out to be at most of order of 
7 g b ~ 10~ 2 T, whereas the thermalization rate of the scalar degrees of freedom receives 
a further suppression C(A 2 /(?|) with respect to 7 g b, A being the quartic self-coupling of 
the Higgs scalar potential and g 2 the SU{2) L gauge coupling. Thus, because of both 
kinematical and loop suppression factors, the thermalization rate of fermions, gauge and 
scalar bosons with the Higgs background turns out to be fairly small. 

Let us now envisage the situation in which the c -background, permeating the thermal 
bath formed by SM degrees of freedom, changes from C = to a non-vanishing value. 
From general point of view the initial field value could be arbitrary with initial distribu- 
tion determined by that. In practise the most of the interesting cases have initially = 
and the applications we consider use that initial value. A realistic example for such a 
situation could be the formation of a subcritical bubble via thermal fluctuations above 
the transition temperature Tf or the passage at a given point of the wall of an expanding 
critical bubble nucleated at Tf in a first-order electroweak phase transition. 

As said above, we can not use the imaginary-time formalism to compute the effective 
Higgs potential for times smaller than ~ 7^ since there is no relation between the density 
matrix of the system and the time evolution operator which is of essential importance 
in the equilibrium case. There is, however, the real-time formalism of Thermo Field 



Dynamics which suites our purposes |13| and is characterized by doubling the degrees 
of freedom of the heat reservoir. It is straightforward to find the (ll)-component of the 
scalar and fermion propagators for the fields 

(13) 

and analogous formulae for the gauge boson propagators. Note, that here /b,f(&) are 
arbitrary distribution functions restricted by requiring the number expectation value to 
be positive. 
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The general expression for the thermal part of the effective potential can be calculated 
by usual way (e.g. by the tadpole method) but the calculation leads now to a different 
expression for the potential than Eq. (]]]) (with / e l q replaced by the generic expression 
/b(f))- The derivation of formula (^]) uses explicitly the concept of thermal equilibrium 
or, as in the case of tadpole method in real time formalism, it is assumed that the 
distribution is a function of energy, i.e. it is <p c dependent. Now, for times smaller than 
7th 1 ) /b,f are wen approximated by 

jeq = (Vlk| 2 +^(T,0 C =O)/T T l)\ (14) 

where m(T, <f> c ) is the plasma mass at finite temperature for each degree of freedom. 

The choice given in Eq. flTjj ) for the distribution functions needs some justification: 
at t < 7^ particles have not had time enough to have real scatterings, and thus the 
distributions have not yet had time enough to feel the change of the background C . So 
they remains the same as they were when the background mass was C = 0. Neverthe- 
less, since forward scattering are much faster than non-forward reaction rates, the free 
dispersion relation of particles get modified and their masses receive plasma corrections 
from the ensemble. Furthermore, since quanta interactions are rather fast, e.g. mediated 
by strong force, the use of the equilibrium distribution function is well motivated. We 
refer the reader to ref. [1141 for more details. 



Calculating the effective potential, we find 

v (t, (pc 1 



2tt 2 



Tr / d\k\ |k| 2 V|k| 2 + m 2 (T> c ) ( e VW^W^)/T T 1 



(15) 



where the trace is taken over all degrees of freedom and m 2 (T, C ) is a general mass 
matrix. Same techniques as used in the calculation of equilibrium potential can be 
applied here resulting a general formula in which the leading terms are included 



V(T, 



2^ 



' m h (T, <p c ) m b (T, 4> c = 0) 



T 



T 



Tr/f 



'm f (T, <p c ) m f (T, C = O)' 



T 



(16) 



The bosonic (fermionic) functions Jb(f) are given by 



h(x, y) 



-7T 



2 1 3 1 

x x 6 H 

3 16 



1 - 7£ + ln(47r) 



4 1 4 i 2 , 

x x lax + 

32 



(17) 



7 



and 

■1 



1 2 1 2 

— 7T y 

24 8 y 



2 1 

x 2 + 
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4 + 7b - ln(vr) 



x 4 + ^-x 4 lnx 2 + ... (18) 



Here 7^ is the Euler constant and m^m is a general mass matrix of bosonic (fermionic) 
degrees of freedom. In the formulas above the dots stand for higher powers of x. It is 
worth to note that in the mass matrices we are able to use the temperature dependent 
masses: the resummation can be done similarly than in the equilibrium case || . This po- 
tential clearly differs from the usual equilibrium potential, whereas the zero-temperature 
part remains naturally unaltered. In particular, the next-to-leading boson contribution 
coming from the infrared region to the cubic term differs from the equilibrum one by a 
factor 2/tt. 

After adding up the two parts we have plotted the effective potential for the Standard 
Model in the particular case of = 65 GeV, Fig. 1. The critical temperature T c 
defined as the temperature at which V(T C , 0) = V(T C , 0+ being the non- vanishing 
minimum, is around 94.27 GeV, higher than the critical temperature for the equilibrium 
case (around 90 GeV) for the same choice of the Higgs mass, see |TJ| for a detail analysis. 

Let us now briefly discuss two physical cases in which our nonequilibrium approach 
to the effective potential might qualitatively effect the conclusions commonly drawn em- 
ploying the equilibrium effective potential. 

As already mentioned, the role played by subcritical bubbles at the onset of a first- 
order electroweak phase transition is still disputed Q. It might be possible that the 
amplitude of thermal fluctuations be so large that the fraction occupied by the symmetric 
minimum in the neighborhood of the critical temperature T c becomes of order unity, thus 
preventing the formation of expanding critical bubbles and hindering any mechanism for 
electroweak baryogenesis. 

To compute the average amplitude of thermal fluctuations it is commonly hypoth- 
esized that the free energy F of a subcritical bubble configuration s b receives contri- 
bution from the equilibrium effective potential. This gives rise to a Boltzmann weight 
~ exp [— F(0 s b)/T] for each configuration. However, subcritical bubbles, being unstable 
objects, tend to shrink. It is rather conceivable that thermal fluctuations are dominated 
by large-size configurations whose lifetime can be roughly estimated to be r s h ~ m~^{T). 
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Since r s h <C 7^ at the critical temperature, particles inside the subcritical bubble do 
not experience thermalization with the s b-configuration. Consequently, the equilibrium 
effective potential is completely inadequate to determine the free energy of such config- 
urations and must be taken over by the nonequilibrium one. 

Our nonequilibrium approach might also have consequences for the determination of 
the velocity v w and width L w of critical bubble walls expanding in the thermal bath at the 
onset of a first-order electroweak phase transition. The determination of these parameters 



has received much attention |L5[ since the discovery of the possibility to generate the 
baryon asymmetry at the electroweak scale. Intuitively speaking, the velocity and the 
shape of the bubble wall depend upon two factors: the pressure difference Ap at the two 
edges of the bubble wall (i.e. between the broken and the unbroken phase), which allows 
for the expansion, and the friction force due to the collisions of the plasma particles off 
the wall. In the previous treatments, the population density /' of each species i has 
been splitted in the equilibrium one plus a small deviation, /' = fL((j) c ) + 8f l . The 
pressure difference Ap is then determined by AV C(l (<j) c ) and the friction force arises due 
to departure from the thermal equilibrium distribution. However, since thermalization 
between the degrees of freedom of the thermal bath and the bubble wall background 
takes a finite time, this picture is correct only if 7^ <C (L w /v w ). If the opposite limit, 
equilibrium is not attained in the vicinity of the bubble wall in the broken phase and the 
nonequilibrium effective potential should play a role in estimating the pressure difference 
Ap between the two phases. 

We are currently investigating both issues [ [Hfl . 
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Figure Captions 

Fig. 1: The Standard Model nonequilibrium effective potential (in units of GeV) for 
the particular choice mu = 65 GeV. 
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